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Self-steepening of optical pulses arises due the dispersive contribution of the (^) Kerr nonlin¬ 
earity. In typical structures this response is on the order of a few femtoseconds with a fixed frequency 
response. In contrast, the effective Kerr nonlinearity in photonic crystal waveguides (PhCWGs) 
is largely determined by the geometrical parameters of the structure and is consequently tunable 
over a wide range. Here we show self-steepening based on group-velocity (group-index) modula¬ 
tion for the first time, giving rise to a new physical mechanism for generating this effect. Further, 
we demonstrate that periodic media such as PhCWGS can exhibit self-steepening coefficients two 
orders of magnitude larger than typical systems. At these magnitudes the self-steepening strongly 
affects the nonlinear pulse dynamics even for picosecond pulses. Due to interaction with additional 
higher-order nonlinearities in the semiconductor materials under consideration, we employ a gener¬ 
alized nonlinear Schrodinger equation numerical model to describe the impact of self-steepening on 
the temporal and spectral properties of the optical pulses in practical systems, including new figures 
of merit. These results provide a theoretical description for recent experimental results presented 
in [Scientific Reports 3, 1100 (2013) and Phys. Rev. A 87, 041802 (2013)]. More generally, these 
observations apply to all periodic media due to the rapid group-velocity variation characteristic of 
these structures. 
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I. INTRODUCTION 

A recent trend in nonlinear optics is the development 
and design of waveguide systems with tunable nonlinear¬ 
ities. In addition to a broad tuning range, these systems 
are characterized by the ability to separate the contri¬ 
butions of the material constituents from the device ge¬ 
ometry This is in contrast to early approaches in both 
glass [H and semiconductors 0 which required a change 
in material composition to modify the waveguide proper¬ 
ties. The role of geometry drastically changed with the 
advent of micro-structured hbers and the demonstration 
that fabrication parameters could be the dominant con¬ 
tribution to the dispersion Q . More recently, it has been 
shown that gas-filled hollow-core fibers can activate or 
suppress nonlinearities such as the Raman effect Q. In 
parallel, rapid advances in integrated semiconductor de¬ 
vices have pushed the forefront of optical science by re¬ 
ducing nonlinear thresholds to sub-femtojoule energy lev¬ 
els [^while simultaneously incorporating dispersion con¬ 
trol p. Among nanostructures, photonic crystal waveg¬ 
uides (PhCWGs) are of extreme interest due to the link 
between geometric fabrication and direct modulation of 
the electric held, giving rise to new physical phenomena 
such as slow-light and enhanced nonlinearity. 

Slow-light refers to light propagating at a reduced 
group-velocity in the medium. Interest in this unique 


property has inspired a large body of research investi¬ 
gating the linear and nonlinear properties of slow-light 
in two-dimensional (2D) PhCWGs over the past decade 
0, [1] ■ Recall the group index rig is related to the waveg¬ 
uide dispersion relation uj{k) and the group-velocity Vg: 
rig = -^ = c^, with frequency w, wavevector k, and 
the speed of light in vacuum c. Of particular signih- 
cance, it was shown that optical effects such as the 
Kerr nonlinearity scale with the group index squared in 
the presence of slow-light Briefly, one factor of rig 
arises from a larger electric field for a given power (non¬ 
linear enhancement), with the second from longer effec¬ 
tive optical path length (linear enhancement). We note 
that the slow-light-enhancement described here is derived 
from the structure. In contrast, material slow light from 
atomic resonances does not exhibit this enhancement Q . 
In PhCWGs we write the effective nonlinear Kerr param¬ 
eter as 7 e// = 7 (^) = (^) ’ ^ith the bulk 

Kerr coefficient n 2 , modal area Ag//, and linear refractive 
index Uq. While the 7 term is well described in the lit¬ 
erature, research into the slow-light enhancement contri¬ 
bution ^ in 2D PhCWGs required significant advances 
in nanofabrication techniques which were only mastered 
the past few years. 

A 2D PhGWG consists of a periodic array of low-index 
dielectric embedded in a high-index material. A com- 



2 


mon experimental configuration which we consider here 
consists of a hexagonal pattern of air holes etched in 
an air-suspended semiconductor slab. Importantly, the 
dispersion of these 2D PhCWGs is highly tunable due 
to selected geometric modifications of the periodic lat¬ 
tice known as dispersion engineering [a, |l0|. The pre¬ 
cise modulation of the waveguide group index enables 
exquisite control over the dispersion and therefore the 
nonlinear properties of the medium. Experimental re¬ 
ports of slow-light enhanced nonlinear Kerr effects in 
2D PhCWGs include demonstrations of solitons, third- 
harmonic generation, and four-wave mixing, amongst 
others [nHil- Despite this strong interest in slow-light, 
the dispersion of the Kerr (w) nonlinearity, or self- 
steepening (SS) term tnl = has received 

surprisingly little attention in these systems. 

The earliest investigations of SS in waveguide systems 
were carried out in glass fiber with these studies em¬ 
phasizing the spectral re-shaping pro per ties Bin or 
the formation of optical shock fronts [l^. Later it was 
shown SS is essential for extending the validity of the 
Generalized Nonlinear Schrodinger Equation (GNLSE) 
envelope approximation down to the single-cycle regime 
[a, and for explaining broadband supercontinuum gen¬ 
eration [l^. In these systems the SS term is determined 
almost exclusively by the wave angular frequency w and 
therefore exhibits a fixed response of about a femtosecond 
for optical frequencies. A more recent numerical inves¬ 
tigation showed the wavelength dependence of the non¬ 
linear Kerr effect in silicon channel waveguides leads to 
values up to tens of femtoseconds near mode cutoff M- 
An alternative approach for generating SS using cascaded 
media as an effective tunable was shown with 
similar strength to traditional media [2^. To date, self- 
steepening in tunable x*'^^ media has not been explicitly 
experimentally demonstrated. 

In this article, we investigate self-steepening in 2D pho¬ 
tonic crystal waveguides. Importantly we show the large 
variation in waveguide group index Ug leads to a new 
physical mechanism for generating self-steepening with a 
characteristic time scale t/vl on the order of hundreds 
of femtoseconds, two orders of magnitude larger than in 
non-periodic waveguide systems [isl; IB • We derive an 
analytic formulation and describe the origin of this ef¬ 
fect. Eurther, we describe structures in which the val¬ 
ues of tjvl are anomalous (negative), hence leading to 
notably different physical effects than previously known 
X^) systems. The broad tuning range of enabled by 
dispersion-engineering make PhGWGs an ideal system 
for further studies of SS. While the magnitude of tnl is 
quite large, the presence of other effects such as group 
velocity dispersion (GVD, P 2 = higher order 

nonlinearities such as multi-photon absorption or free- 
carrier effects can disrupt the ideal dynamics. We conse¬ 
quently describe the experimental situations in which SS 
is expected to contribute significantly in the semiconduc¬ 
tor system under consideration using a numerical model. 
This analysis supports recent experimental results show¬ 


ing pulse temporal advance in PhCWGs [ITI,!!^. Though 
the Raman-effect is narrowband in semiconductors and 
negligible here, this analysis could be extended to in¬ 
clude periodic glass media where Raman and Brillouin 
effects must be considered [B- These results provide 
the first theoretical description of giant and tunable self- 
steepening in nanoscale optical waveguides. 

More generally, this investigation applies to all periodic 
media (ID, 2D, 3D) where a dispersive arises 

due to strong group-index modulation near the band- 
edge. Given the importance of the dispersive nonlinearity 
tjvl in explaining supercontinuum broadening in fibers, 
we expect the new terms elucidated here to be critical 
for accurately describing this phenomenon in photonic 
crystals. 


II. SELF-STEEPENING IN PHCWGS 
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FIG. 1. Group index curves and self-steepening parameters 
of PhC waveguides, (a) We show three typical waveguides: 
(b,green) a dispersion-engineered structure with a quasi flat 
plateau, (c,cyan) a W1 line-defect waveguide, and (d,blue) 
and a dispersion-engineered waveguide with a peak, (b)-(d) 
Show the self-steepening parameter tnl = 71 / 7 e// (black 
line) as a function of wavelength for the three structures. The 
different contributions to SS are shown: l/w(black-dashed), 
Sdu;rig/ug (red) and duiAf.ff /Af.ff {blue). 
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Though near-arbit rary dispersion profiles are possible 
in periodic media focus on three specific 

experimentally demonstrated structures for clarity. Fig¬ 
ured] shows three group index curves for PhCWGs with 
different dispersion relations: (i) a standard line defect 
waveguide of one missing row of holes in a hexagonal 
lattice (Wl), (ii) a dispersion engineered waveguide @ 
exhibiting a plateau, and (iii) a dispersion with a pro¬ 
nounced group index peak. The dot indicates a point of 
interest we investigate in this work. 

If we assume the Kerr nonlinearity evolves linearly then 
the impact of dispersive nonlinearity on pulse propaga¬ 
tion dynamics is modeled by adding a first order Kerr 
correction term to the GNLSE [l3, [l^ : 

dzA + ^P2dttA = f7e//(l + iTNLdt)\AfA. (1) 

Here t is the relative time in the reference frame of the 
pulse with A{z,t) = y/Po{z, the electric-field 

envelope with power Po and phase </>. The last term 
is referred to the self-steepening or shock term and ex¬ 
pressed as: 


contributions on this scale. Second, the sign of tjvl is 
negative. Gonsequently the spectral and temporal prop¬ 
erties of the nonlinear waves behave in an opposite or 
anomalous manner as we will examine. Third, the dom¬ 
inant contribution arises from the dispersion term (red) 
with a modification in the opposite direction due to Agf j 
(blue). If we ignore the area contribution, we approxi¬ 
mate TnL ^ _ 2 _ 9 p_ _ 2 c§ 2 ^ 

^ Ug OUJ rig 

Examining the characteristics of the waveguides indi¬ 
vidually, we hnd different trends for each. Regarding 
the Wl waveguide Fig[TKc), steadily decreases ap¬ 
proaching the mode cutoff (increasing wavelength). A 
value of about -200 fs is obtained close the the band 
edge, however in a region where propagation loss is large 
in this type of PhCWG [1^. On the contrary for the 
dispersion-engineered waveguides the values are mostly 
negative with a small positive region for the structures 
presented. Notice in this case values as large as -200 
fs (Fig HKb)) and -400 fs (Fig [I^d)) are reached away 
from the band edge where propagation losses are reduced 
[^tL [ 2 ^ . The lower linear loss of the latter structures has 
implications for practical observation of these effects. 


^ ^ _ 1 dAgff 2 dug 

ygjy w n2 dui Agff duj Ug dco ’ 

where we used the relationship 71 = d^'yeff. The two 
first terms exist in bulk material and contribute up to 
a few femtoseconds. In practice the Kerr dispersion is 
often neglected. These are the traditional self-steepening 
terms known in unstructured waveguides and do not play 
a role here with these small magnitudes M- The third 
term encompassing effective modal areas Agff was con¬ 
sidered in fibers [^, and theoretically in silicon channel 
waveguides with surprisingly little attention in PCFs 
[ 2 ^. Further, in those earlier works the area term con¬ 
tributed a few percent whereas here the slow-light modes 
exhibit rapid variation in spatial profile with ui and this 
term contributes approximately ~ 25 % to the total tnl- 

The key physical insight of this work is the realization 
that the final term due to the dispersive group-index Ug 
gives rise to a previously unknown mechanism for control¬ 
ling self-steepening. In contrast to earlier observations, 
in our photonic crystal waveguide the magnitude of tnl 
is set by a combination of the group-index Ug, GVD (/I 2 ) 
and effective modal area Agff. As these three parame¬ 
ters depend strongly on the geometry, we clearly see the 
advantage of using nanostructures to study the dispersive 
nonlinearity. We now examine this effect in detail. 

Figures [TKb)-(d) show the evolution of the ratio t^l 
for the waveguides in Fig[TJa). The contribution of the 
different terms composing Eqn. [5] are also represented. 
We note three major observations unique to the PhGWG 
system. First, is two orders of magnitude larger 
than unstructured waveguides (black dashed line) where 
tnl = ^ ~ 1 fe. Notice the ^ contribution appears to be 
near-zero and completely flat compared to the PhGWG 


III. TEMPORAL AND SPECTRAL 
PROPERTIES DUE TO ANOMALOUS 
SELF-STEEPENING IN PHCWGS 

We now describe the physical implications of the 
self-steepening term on nonlinear wave propagation in 
PhCWGs. For that purpose we consider typic al p arame- 
ters found in recent nonlinear experiments [Tl|, [U, [ 2 ^, 

[ 3 ^. We take 7 e// = 1600 (W.m)“^ (ng=15, no=3.17 for 
GalnP), an anomalous dispersion of /32 = -7.7 ps^/mm, 
^ = 6xlO“^®m^/W and modal area Agff = 0.34 
[g. The dispersive nonlinearity is tml= -220 fs as 
detailed above. The input pulses are TpwHM = 2.3 
ps (full-width at half-maximum of a hyperbolic secant, 
Tfwhm = 1-76 To) with Pq = 3 -10 W (6 - 20 pJ/pulse). 

rp‘2 

The dispersion length Lp = -^ is computed as 220 fim. 
Importantly, throughout this work we purposely main¬ 
tain small soliton numbers (N < 2) so as to avoid more 
complicated soliton dynamics modulating the peak in¬ 
tensity and pulse duration (2ll |. 

At this point we are focusing on the basic physical 
effects resulting from the unique photonic crystal disper¬ 
sion in the ‘ideal’ system. In the next section we will in¬ 
troduce the full effects present in typical semiconductor 
waveguides and describe how these results are modified. 
While based on actual experimental structures, note that 
the conditions may not be optimal for emphasizing the 
self-steepening effect and we invite the community to ex¬ 
plore the parameter space further. 

The normalized self-steepening parameter s is: 


For our parameters s= -0.1, more than five times larger 
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than non-periodic waveguides [l5l | . This is even more re¬ 
markable when one considers the pulses are 2.3 ps long 
compared to the sub-100 fs pulses required in unstruc¬ 
tured media where s = The large value of s re- 

UllQ ^ 

quires much shorter length scales to observe the asso¬ 
ciated effects of self-steepening. A pulse experiencing 
self-steepening will eventually develop a shock front af¬ 
ter propagating a shock lenathfld^ of about: 

z, = 0.43^, (4) 

|s| 

where L]\il = {leffPo)~^ is the nonlinear effective 
length, and the numerical constant depends on the ac¬ 
tual pulse shape, here a hyperbolic secant [l^. Typical 
shock lengths in our structures are Zg sa 350 ^m for a 2.3 
ps pulse with peak power of 8 W. As this article is fo¬ 
cused on self-steepening, we do not describe the physics 
of optical shock waves in detail here. Nonetheless, this is 
a well known and useful length scale for estimating the 
relative scaling of self-steepening which we adopt here for 
convenience. Note that the mechanism presented here is 
not the only possibility for developing shock fronts. A 
highly nonlinear medium in presence of weak normal dis¬ 
persion, for example, could also lead to shock formation 
even though no dispersive nonlinearity is present (3lj |. 

The nonlinear dispersion tnl plays a key role in the 
pulse dynamics of both the temporal and spectral prop¬ 
erties. We first address the impact of SS on temporal 
shape and delay. A subtle point that must be addressed 
is the dual role of /32- First, it has been shown in earlier 
work that P 2 dissipates shock fronts [T^. Second, and 
separate to this point, we showed above that P 2 is intrin¬ 
sically linked to the large magnitude of tjvl. As a result, 
one cannot ignore f32 for a SS effect arising from a strong 
modulation of the group index Ug and consequently ob¬ 
serving a tilted line-shape is unlikely for SS derived from 
this method. While the temporal shape is not modulated 
in this case, the temporal arrival time is affected as we 
show below. 

Figure HJa) shows the temporal profile obtained for 
the point indicated by the blue dot in Fig [Ua-c) after a 
propagation distance of 200 fim. The input pulse (black 
line) has a peak power of 8 W {Lnl=80 /im), hence 
the pulse has completed about two and a half nonlinear 
lengths Lnl- The thick-red line shows the case where 
we include only the Kerr and SS contributions to Eqn. [T] 
That is, we include the /?2 contribution to tnl but neglect 
the temporal dispersion term dtt- Notice that since s 
is negative here the pulse peak tilts forwards in time, 
which is opposite to earlier studies with a positive self- 
steepening term. Moreover, with the large s value here, 
a steep temporal leading edge is already clearly visible 
after a propagation of a just few Lnl- Simulations with 
all the terms in Eqn. [1] are shown as the dashed blue 
line. In contrast to the ideal dispersion-less case (thick 
red) where the temporal trace exhibits a clear abrupt 
leading edge characteristic of shock formation, the shock 
front is less pronounced for the same propagation length 


and the pulse tends to preserve its initial shape (soliton 
effect). 

A separate temporal effect in the presence of self- 
steepening is a shift in pulse arrival time. Assuming the 
pulse duration Tq > and a moderate soliton number, 
the dispersive nonlinearity acts as a small perturbation 
that slightly modifies the group velocity AT = j^tnl 
[ll,|3l,[33- Here with the negative tnl value, the pulses 
are expected to advance in time, once again in contrast 
to the delay of earlier observations. This effect is clearly 
visible by the temporal advance of the two traces exam¬ 
ined in Eig [21 One of the most challenging aspects of 
observing self-steepening in real PhCWGs is the strong 
resemblance of SS and free-carrier dispersion (ECD) in 
the time domain. We investigate this in detail below. 

Eigure|2Kb) shows typical pulse spectra for the negative 
s values characteristic of PhCWGs. The pulse spectrum 
is clearly rendered asymmetric in the presence of self- 
steepening. The blue-shifted peaks become more intense, 
while the red components are notably broader compared 
to the symmetric broadening characteristic of SPM. The 
role of P 2 is less pronounced. 
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(b) 
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FIG. 2. Temporal and spectral pulse properties of anoma¬ 
lous self-steepening in PhC waveguides. (a) Temporal 
pulse shapes, (b) Spectral pulse shapes. Notice the be¬ 
havior is opposite to unstructured media exhibiting nor¬ 
mal self-steepening. Traces are shown for propagation of 
200 fim « 2.5Lnl- Legend: Kerr and SS only (thick red), 
full Eqn. [T] (dashed blue), and input pulse (black). 


IV. SELF-STEEPENING IN 
SEMICONDUCTORS 

In our analysis thus far we have ignored several ef¬ 
fects present in practical systems. The role of linear 
loss has been previously studied and in the limit azg >I 
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shock waves are unobservable altogether [I^. For the 
waveguide experiments referenced above the linear loss 
is given by a = 15 dB/cm (at ng = 15). This yields 
azs = 0 . 11 , indicating these waveguides can in principle 
support shocks. 

However, in semiconductor media we must also con¬ 
sider nonlinear absorption of multiple photons across the 
electronic bandgap. For a typical wavelength of 1550 nm 
(photon energy of 0.8 eV), silicon {Eg = 1.1 eV) is re¬ 
stricted by two-photon absorption (TPA) and the wide- 
gap material GalnP {Eg =1.9 eV) is limited by three- 
photon absorption (SPA). While the optical properties 
of silicon have been widely studied, only over the past 
years we have investigated the properties of GalnP 
and established this material as a viable platform for non¬ 
linear optics at 1.5 /rm [n|, [2l|,[s^,[13]■ In simplest terms, 
nonlinear absorption damp the dynamics similar to lin¬ 
ear absorption. An order of magnitude estimate shows 
loss due to TPA (SPA) requires a 2 lZs < 1 {a^Pzg < 1), 
with the intensity I = , to observe a shock. Physi¬ 

cally these ratios compare the strength of the nonlinear 
loss to the self-steepening term. 

For the TPA case (e.g. Silicon, 02 = 1 GW/cm (S5j |l 
the ratio a 2 lzs = ss 1.8 indicating SS-induced 

shocks are not generally accessible in this material, even 
for this large s= -0.1. Here we have defined a new nonlin¬ 
ear figure-of-merit (FOM) for self-steepening and TPA. 
Its form is noticeably similar to the well-known version 
for Kerr-TPA switching {-^^) with the additional 
term for SS. Notice this ratio is independent of power 
for TPA and therefore SS will always be much weaker 
than TPA. Since TPA dominates the SS effect, we focus 
mainly on the SPA system in the following analysis. 

The SPA material (e.g. GalnP, a3=0.01S GW^/cm^ 
[S7| l is much more favorable, yielding a^l'^Zg = 
0 ^ ^ 2 which is satisfied for intensities up to about 

|s| koUQ, — 

~50 GW/cm^. This threshold is intensity-dependent due 
to the different nonlinear orders of Kerr and 
SPA. Including the slow-light enhancements 03 oc and 

712 C)C ni this relation becomes / 11:9 ) < ^ and is 

9 |s| kon2 \n.o J — 

satisfied up to «10 GW/cm^ for the conditions here. In 
contrast, the TPA-limited SS case does not scale with 
slow-light. Note these estimates are only indicative and 
do not correspond to strict thresholds. For larger pow¬ 
ers, free-carrier absorption would play an important role 
until the peak power falls below its threshold. Now that 
we have established the role of nonlinear loss, we describe 
the free-carrier effects. 

Free-carriers generated via these nonlinear absorption 
mechanisms have an equally significant impact in the 
pulse dynamics through both dispersive (FCD, npc) and 
absorptive (FGA, a) contributions. Importantly, the 
physical signatures of anomalous self-steepening strongly 
resemble those of FCD, especially in the time domain. 
Namely, like anomalous SS, the pulse also advances as 
a function of input power due to FCD combined with 
anomalous GVD as our recent experiments show [13, . 


A broader GNLSE including all of these effects is: 


dA a ,^2 52 A 


-I- {ikoUFC - 
-^\A\^A. (5) 


For the modeling below, we take the slow-light scaled 
values of the bulk parameters as described in prior 
literature [ 13 , [13 ■ The parameters are a^eff = 

(^)'= 18 m-bW-^ npc,eff=npc (^)= -1.8 
Xl0-26m3, cTe// = 1-3 xlQ-^W. The 3PA 

area is A5e//[13|- We have also included the third- 

order dispersion (TOD, P 3 , = = -1-0.7 ps^/mm) 

for completeness as this cannot be ignored in the real 
system. Regarding one case briefly presented below 
where the material is limited by TPA (silicon), we take 

a 2 eff = (ttY) ~ (W.m)“^. It is worthwhile 

to point out that these parameters correspond to actual 
experimental parameters and therefore are immediately 




FIG. 3. Temporal and spectral behavior of anomalous self- 
steepening in after a propagation distance of 500 pm in realis¬ 
tic PhCWGs. (a) Temporal traces for 3PA-limited materials: 
full (green), FGD alone (blue), self-steepening alone (red). 
The black lines correspond to the input pulse, (b) Spectra 
for 3PA-limited materials, (c)-(d) Pulse temporal advance as 
a function of power for (c) 3PA-limited (solid) and (d) TPA- 
limited (dashed) systems. Green curves correspond to the full 
model in Eqn. [S] Blue lines are the case where we neglect SS, 
whereas red lines indicate when FCD is neglected. We include 
the full 3PA result in (d) to compare the relative scales of the 
two systems. 
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realizable in current systems. 

Figure [31Ja) shows the pulse temporal shift due to the 
competing self-steepening and FCD effects in the SPA 
limited material (a 2 = 0) at a peak power of 8 W. When 
considered in isolation the FCD-GVD curve (blue, no SS) 
and SS-only curve (red, no FCD) each would contribute 
a few picoseconds of delay. Moreover they have a rela¬ 
tively similar temporal shape. Note these effects do not 
add linearly, but rather compete for power and interact 
dynamically to yield the full result (green). 

Importantly, the spectral features of FCD and SS are 
distinct. Figure |3Kb) shows the spectral properties of 
pulses propagating in the SPA system. Considering only 
the self-steepening effect (red curve) results in a minor 
modulation to the symmetric shape expected from SPM- 
only and does not shift the spectral center-of-mass. In 
contrast, the FCD induces a clear blue shift in the pulse 
center-of-mass [2l|. When we consider these effects si¬ 
multaneously (green curve), the FCD is clearly the dom¬ 
inant contribution. Thus these effects are more easily 
distinguished in the spectral domain. 

The nonlinear scaling laws of SS and the FCD- 
CVD temporal shift are also notably different. No¬ 
tice tk!l oc Po whereas FCD^^^ oc (TPA) or 

FCD^^^ oc Pq (3PA), where we have written FCD^™), 
with m indicating the order of nonlinear absorption gen¬ 
erating the free-carriers. Figures EKc) and (d) report the 
scaling of SS and FCD as a function of power for (c) 
SPA and (d) TPA-limited materials. We observe SS (red 
curve) scales much more slowly compared to FCD (blue 
curve). The full simulation (green dashed) more closely 
follows SS at low powers and the FCD trend at higher 
power. For the TPA case shown in (d) the SS induced de¬ 
lay is noticeably smaller due to stronger nonlinear TPA 
loss which caps the peak power more than the SPA case. 
In (d) we also show the SPA curve for comparison, high¬ 
lighting the greater temporal shift in this case. 

The striking similarity of the temporal advance and 
pulse shapes from both SS and FCD-CVD are highlighted 
in recent experimental results. The CalnP waveguides 
in both cases are similar to that in Figjljc), with rela¬ 
tively small values of t^l compared to that highlighted 
in this work. In a ThPA system the temporal advances 
are attributed to FCD-CVD supported with numerical 
modelling , however with no SS term. Similar results 
were shown in Refs.dl,!!^ in silicon. In contrast, a pulse 
advance attributed to SS only was reported in Ref. . 
Critically that report did not include the important con¬ 
tribution of FCD-CVD term, but rather attributed it to 
SS alone. As we have shown here, the FCD-CVD ad¬ 
vance plays an equally important role as SS and cannot 
be ignored. 

The true physical situation is likely a combination of 
these effects though the exact scaling would be challeng¬ 
ing due to soliton dynamics modulating the peak power 
and requires careful experiments and numerical modeling 
to discern the effects. Nonetheless, one could use the mo¬ 
ments method in Ref. to estimate the magnitude of 


the temporal shift from FCD-CVD or the SS time shift 
equation given above for order-of-magnitude values. 

One last important aspect to consider is the relative 
impact of the instantaneous SS and non-local FCD ef¬ 
fects on pulse delay during propagation. The SS induced 
delay depends directly on the instantaneous peak power 
and hence will eventually decrease due to linear and non¬ 
linear loss. In contrast, the effect of FCD on the group 
velocity persists even if the pulse peak power decreases. 
Physically this results from the spectral FCD frequency 
blue-shift and accompanying change in the group-velocity 
experienced by the pulse |4Cll |. 

Finally we show the evolution of the nonlinear dynam¬ 
ics of the pulse propagating down the waveguide. Figure 
SKa) shows the change in pulse delay along the waveg¬ 
uide length for the two mechanisms in the SPA system 
at a peak power of 8 W. The change in delay due to the 
SS effect (red, FCD=0 and FCA=0) is strongest in the 
first 200 /im and then tapers off due to loss. The FCD 
effect (blue, 71 = 0 ) continues to experience a change in 
delay right until the end of the propagation with the full 
curve (green) more closely resembling this case. Figure 
SKb) shows the corresponding peak power evolution which 
gives insight into the dominant nonlinear mechanism as 
the pulse evolves. We attribute the initial increase in 
peak power to soliton compression. 
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FIG. 4. Role of instantaneous SS and non-local FCD on 
pulse delay, (a) Contribution of SS alone (red) and FCD 
alone (blue) to the total pulse delay (green) for propagation 
lengthts up to 1 mm. (b) Peak power evolution correspond¬ 
ing to the cases presented in the upper panel. Note the low 
overall peak power after 500 fim of propagation. Legend: full 
model (green), no free carrier effects (red), no SS (blue). 

V. CONCLUSION 

In this Letter, we investigated the nonlinear self- 
steepening effect in photonic crystal waveguides. Our 
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first principles derivation in the nanostructured periodic 
waveguides revealed a self-steepening term two orders 
of magnitude larger than typical systems. Importantly 
the self steepening coefficient is determined by 

the geometric parameters of the waveguide, offering a 
large tuning range of both positive and negative values. 
The origin of this giant tjvl is the strong dispersion of 
PhCWGs counterbalanced by a modal area contribntion. 
We considered the role of higher-order effects in practical 
systems with new figures of merit, concluding that the 
nonlinear loss quenches the dynamics. We showed that 
the principal physical signature of the anomalous SS 
effect is a temporal forward tilt and pulse advance in 
contrast to the delay observed in normal SS media. In 


the semiconductor waveguides these effects compete 
with FCD which also advance the pulse, whereas the 
spectral signatures of self-steepening and FCD are 
distinct. We suggest future experiments exploring the 
full range of these dynamics be undertaken to reveal 
the full dynamics of this giant tunable self-steepening 
mechanism, especially in the supercontinuum regime. 
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